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Abstract 

Using mixed analytical and numerical methods we investigate the development of singularities in 
the heat flow for corotational harmonic maps from the d-dimensional sphere to itself for 3 < d < 6. 
By gluing together shrinking and expanding asymptotically self-similar solutions we construct 
global weak solutions which are smooth everywhere except for a sequence of times T\ < T2 < 
■ ■ ■ < Tf. < 00 at which there occurs the type I blow-up at one of the poles of the sphere. We give 
evidence that in the generic case the continuation beyond blow-up is unique, the topological degree 
of the map changes by one at each blow-up time Tj, and eventually the solution comes to rest at 
the zero energy constant map. 
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INTRODUCTION 



Let M and N be Riemannian manifolds with metric tensors and Gab m some local 
coordinates {x 1 } and {X" 4 }. A map X : M — > N is called harmonic if it is a critical point 
of the energy 

/i f)x A f)X B 

In this paper we consider harmonic maps from the d- dimensional unit sphere to itself, i.e. 
M = N = S d with gij and Gab being standard round metrics. We parametrize S d by 
spherical coordinates (9, 0), where 9 is colatitude (0 < 9 < it) and is a point on the equator 
S^ 1 of S d . We restrict our attention to corotational maps of the form (9,<p) — > (U(0),<p). 
For such maps we have 

where for convenience we dropped the multiplicative factor v o/(5' d_1 ) coming from the inte- 
gration over 0. The Euler-Lagrange equation corresponding to the energy (j2J) reads 
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It was shown in []J that for 3 < d < 6 Eq.(j3]) has a countable sequence {U n } of smooth 
solutions of degree zero and one. These solutions may be viewed as excitations of the ground 
states: the constant map Uq — and the identity map Lq = 0, for even and odd values of 
n, respectively (for d > 7 these excitations disappear). For n — > oo the solutions U n {9) 
converge (nonuniformly) to the (singular) equator map = n/2. Later, Corlette and 
Wald [2| rederived and extended these results using Morse theory methods. Their approach 
helped to identify the two key features which are responsible for the existence of infinitely 
many solutions: the presence of the antipodal reflection symmetry U — > n — U and the 
existence of the singular map = 7r/2 of infinite index which is invariant under this 
symmetry. An essential ingredient of the Morse theoretic argument is an energy decreasing 
flow in the space of maps. In Q] this flow was defined in a somewhat ad hoc manner to 
ensure that it has all the desired technical properties. One might wonder if it is possible to 
repeat the Corlette- Wald argument using the ordinary heat flow. This would be interesting, 
for instance, in numerical implementations of the argument for similar systems. The main 
technical difficulty is that the heat flow can develop singularities in finite time. If this 



happens, in order to save the argument, one must find a way to continue the flow past a 
singularity in a unique manner. Although an analysis of this issue was the original motivation 
for this paper, the problem of uniqueness of continuation beyond blow-up in the heat flow 
for harmonic maps seems interesting in its own right, regardless of possible applications to 
elliptic problems. 

The aim of this paper is two- fold. First, we describe the precise asymptotics of blow-up 
in the heat flow for corotational harmonic maps from the d- dimensional sphere to itself for 
3 < d < 6. We show that blow-up has the form of a shrinking self-similar solution (shrinker, 
for short). In turns out that among infinitely many shrinkers (whose existence was proved 
by Fan [3|), there is exactly one which is linearly stable. We provide numerical evidence 
that this stable shrinker determines the generic profile of blow-up. Second, we continue the 
flow past the singularity by gluing a suitable expanding self-similar solution (expander, for 
short). We find that there is exactly one expander which can be glued to the stable shrinker 
and consequently the continuation beyond the generic blow-up is unique. 

The scenario of incomplete blow-up and self-similar global "peaking solutions", that is 
solutions which shrink self-similar ly, blow up, and then expand self-similar ly for a while 
(with this scenario possibly repeating a number of times) has been studied in the past 
: or the harmonic map flow |4J and other parabolic equations: the semilinear heat equations 
5, 6], the mean curvature flow {4, 7-9], the Yang- Mills flow 10], the Ricci flow 11], and more 
recently for a fourth-order react ion- diffusion equation [l2|. Most of these studies emphasized 
non-uniqueness of continuation beyond blow-up. To our knowledge, this is the first work 
which demonstrates (by heuristic and numerical means) that for the generic blow-up the 
continuation is unique. As we shall see below, the uniqueness of continuation is contingent 
upon certain quantitative properties of self-similar solutions and thus may be hard to prove. 

The rest of the paper is organized as follows. In section 2 we introduce the heat flow for 

equivariant harmonic maps from S d (or M. d ) into S d and recall basics facts about blow-up. 

Section 3 is devoted to self-similar solutions of the heat flow for harmonic maps from M. d to 

S . Using matched asymptotics we derive asymptotic scaling formulae for the parameters 

of self-similar solutions. In section 4 we analyze the linear stability of self-similar solutions. 

In section 5 we study the continuation beyond blow-up and formulate the main result of 

this paper, that is the conjecture about the uniqueness of continuation in the generic case. 

Numerical evidence supporting this conjecture is presented in section 6. Finally, in section 7 

we indicate possible extensions of our results. 
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PRELIMINARIES 



We consider the heat flow equation 
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with initial and boundary conditions 

U (0, 9) = h(9) G C°°[0, 7r] , (5) 
U(t, 0) = h(0) = , (6) 
E7(t,7r) = /i(7r)=ibr, (7) 

where an integer is the topological degree of the map. As long as the flow is smooth, the 
solution remains in the given homotopy class (i.e., the degree k does not change). It follows 
from (HI) that for a smooth solution there holds 



dE 
~db 



Iff sin^e d6 , (8) 



which shows that Eq. (J3J) is the gradient flow for the energy (TJ]) . Thus, one might expect that 
for t — > oo the solution U (t, 9) will converge to a critical point of E, i.e. a harmonic map. 
Unfortunately, as mentioned in the introduction, this expectation is too naive because in 
general the flow develops singularities in finite time. Indeed, it follows from general results 
for harmonic maps between compact manifolds that for any initial map with nonzero degree 



and sufficiently small energy the solution must blow up in finite time (see Thm 1.12 in 131]) . 

For equation (J3J), by symmetry, the singularity must occur at one of the poles. Since the 
blow-up is a localized phenomenon, the curvature of the domain manifold plays no role in 
the description of asymptotics of blow-up. Thus, from here until section 6 we replace the 
domain S d by its tangent space at the pole, M d , and consider the heat flow for corotational 
harmonic maps from R d to S d 

^ = ^i(^"V) r -^sin(2 M ), (9) 

where u = u(t,r) (r = \x\). Such maps enjoy scale invariance: if u(t,r) is a solution, so 
is u\(t,r) = u(t/X 2 ,r/X) for any positive number A. Solutions which are invariant under 
rescaling, that is u\ = u, are called self-similar. The self-similar solutions play the key role 



in the dynamics of type I blow-upQ so the next three sections are devoted to their existence 
and properties. 

Throughout the rest of this paper we assume that 3 < d < 6. 

SELF-SIMILAR SOLUTIONS 
Shrinkers 

Let us assume that a solution of Eq. ([9]) develops a type I singularity at r = in a finite time 
T, i.e. (T — t)uf,(t, 0) is bounded as t /* T. To describe the formation of the singularity it 
is convenient to introduce new variables 

s= _l n ( T _ t ) ? y= _Z_ ) f(s,y) = u(t,r). (10) 
In these variables Eq.flH]) takes the form 

fs = \ (Pfy) y ~ ^ sin(2/) , p(y) = y d ~ l exp(-y 2 /4) . (11) 

p z y 

This equation can be viewed as the gradient flow for the functional 



z Jo \ y 

We shall refer to S(f) as the conformal energy because it is the energy for maps from 
(M. d , e 2 ( d - 2 ) 5) to S . A simple calculation gives 

ds r f2 



f/s , ft pdy, (13) 

hence the conformal energy is monotonically decreasing. The assumption that the blow-up 
is of type I implies that f y is uniformly bounded as s — > oo, hence the flow must converge 
to a critical point of the conformal energy, that is a solution of the Euler-Lagrange equation 
8£{f) = for f(y) 



/'M^-f)/'-^(2/H0. (14) 



Note that an endpoint of evolution cannot be the trivial solution / = as this would 
contradict the occurrence of blow-up at time T. Thus, the study of type I blow-up reduces 
to the study of nonconstant solutions of Eq. (fl4|) . We shall call such solutions shrinkers. 



1 It is customary to divide singularities into two types: a singularity for which (T — i)|Vu| 2 is bounded as 
t /*• T is said to be of type I; otherwise it is said to be of type II. 
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Let us discuss now existence and properties of shrinkers. Regular solutions of Eq.flT 
behave near y = as follows 



f{y) = ay 



a{Ada 2 - Aa 2 - 3) 



y 3 + 0{y z 



12(2 + d) 

where a is a free parameter. Regular solutions at infinity behave as 

(d - 1) sin(26) 



(15) 



f(y) 



7T 

2 



b- 



2y 2 



+ 0(y- 



(16) 



where b is a free parameter. Using a shooting method Fan 3] proved that for 3 < d < 6 
there is an infinite sequence of pairs (a n , b n ) for which the local solutions (jl~5"I) and (Fl6|) 
are smoothly connected by a globally regular solution f n {y)- The integer index n denotes 
the number of intersections of the solution f n (y) with 7r/2 (see Figure 1). As n — > oo 
the shrinkers converge (nonuniformly) to the equator map /oo = tt/2 and correspondingly 
E(f n ) — > E^foo) = 2 d ~ 1 r(^-). Some quantitative characteristics of shrinkers are displayed 
in Table I. 



n 


a n 


b n 




1 


2.738753 


0.573141 


1.485688 


2 


2.927644 • 10 1 


-0.184519 


1.738165 


3 


3.141830 • 10 2 


0.566142 • 10" 1 


1.771588 


4 


3.376630 • 10 3 


-0.172776 • lO^ 1 


1.776470 


5 


3.629513 • 10 4 


0.527011 • 10~ 2 


1.778116 


6 


3.901390 • 10 5 


-0.160744 • 10~ 2 


1.779706 


7 


4.193637 • 10 6 


0.490287 • 10~ 3 


1.781650 


8 


4.507777 • 10 7 


-0.149542 • 10" 3 


1.784095 


9 


4.845449 • 10 8 


0.456120 • 10~ 4 


1.787199 


10 


5.208415 • 10 9 


-0.139121 • 10" 4 


1.791128 



TABLE I. Parameters of the first ten shrinkers for d = 3. 



From the shooting argument in [3j it follows that a n — > oo and b n — > as n — > oo. We 
shall now use this fact to describe the behaviour of shrinkers for large n (we follow here a 
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FIG. 1. The profiles of the first four shrinkers for d = 3. 
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similar argument given in Let £ = ay and 0(£) = f(y). In terms of these variables 
Eq. ffMl) becomes 



with the initial condition 0(£) ~ £ near £ = 0. For a — > oo, solutions of this equation tend 
uniformly on any compact interval to solutions of the limiting equation 

0" + ^0'-^sin(20) = O (18) 

with the same initial condition ~ £ near £ = 0. Using the standard phase-plane analysis 
we get for 1 -C £ <C a 



7/^ 71 * d -' 2 . , i „ rx V8rf — ^ — 8 

0(0 - 2 +«r~sm(wln£ + 5), w= , (19) 

where the amplitude a and the phase 6 are uniquely determined by the initial condition 
0'(O) = 1. Returning to the original variables we obtain for 1/a <C ?/ < 1 

f(y) ~ — + ay"^ sin(w In y + u In a + 5) . (20) 

On the other hand, for y ^> 1/a the solution is close to 7r/2 so we can write 

f(y)~l + bh(y), (21) 
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where h(y) is the solution of the linearized equation 



— h!+*=±h = 



V 2y y 
normalized by the condition h(oo) = 1. For 1/a < t/ C 1 we have 



(22) 



h(y) aiy 2 sin (a; In y + 5i) 



(23) 



where «i and 5i are uniquely determined by the normalization condition h(oo) = 1. Using 
and matching the solutions (|2"U1) and f|2TT) we get 



a 2 a sin (a; In ?/ + win a + 5) ~ feai sin(a;lny + 5^ 



(24) 



hence 



a; In a + 5 ~ Si + nn, 
which yields the scaling laws for large n 

a n ~Cexp(^, b n ~ (-l)"L>exp 



a 



6 ~ (-1)"— a"" 

Q!i 



n(d-2)n 
2ui 



(25) 



(26) 



where C = exp[(<5i — S)/u\ and D = — C d ^ . Numerical parameters of shrinkers are 
displayed in Table 1 and compared with the asymptotic expressions ( 1261) in Fig. 2. 
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FIG. 2. The asymptotic formulae (I26p are shown to give excellent approximations for the param- 
eters of shrinkers even for small n (here d = 3). 
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Expanders 



To describe the behaviour of solutions for t > T we introduce new variables 

a = \n(t-T), y = -L= } F(a,y) = u(t,r), (27) 
in which Eq.([9]) takes the form 

Fa = ^ (RF y ) y - ^ sin(2F) , R(y) = y^ 1 exp(y 2 /4) . (28) 
We shall refer to time-independent solutions of this equation as expanders. Eg. (|28p have 



been very recently studied by Germain and Rupflin [14j who established interesting results 
concerning existence, uniqueness, and stability of expanders. Below we complement these 
results by a more detailed formal quantitative analysis (which is essential for our purposes). 
Expanders satisfy the ordinary differential equation 



i ,» + (l_I + |) F '_l_l sin(2F) = (29) 



with the regularity condition F(y) ~ Ay near y — 0, where A is a free parameter. In 
contrast to shrinkers, expanders are globally regular for any A. This is due to the strong 
damping term |F' in ( 129]) which drives F'(y) rapidly to zero as y — > oo and guarantees that 
lim^oo F{y) exists. Let B = lim^oo F(y) — n/2. It is routine to show that B depends 
continuously on A. In order to get a more precise asymptotic behaviour, we rewrite Eq.( l29l 
in the integral form 



F>(y) = ^—y 1 - d e - y2/A S d -V 2/4 sin(2F(s))cfe, (30) 



o 



d-1 1 ■ 

yy) = 

and compute the limit 

* „ s , , fr? s d ~ 3 e s2 / 4 sm(2F(s))ds , 
lim y 3 F'(y) = (d - 1) lim ^ d _ A 1 = - (d - 1) sin(25) , (31) 

where the last step follows from rHopital's rule. Therefore, if B ^ 0, we have for large y 

F(y) = l + B + ^r Mm + 0(y-*) . (32) 

We note in passing that Eq.f PBl) below implies that there is an infinite countable subset of 
parameter values for which B(A) = and 

F(y) -tt/2 ~ cy~ d e- y2/A as y ->■ oo . (33) 



The variational proof (using a renormalized energy) of existence of such rapidly decaying 



expanders was recently given in 14( . Since these solutions do not seem to participate in the 
dynamics of blow-up, we do not pursue them here in more detail. 

Next, we derive asymptotic approximations of the function B(A) for small and large 
arguments. For small A we substitute F(y) = AF(y) into Eq. fl29l) and take the limit A — > 
to obtain the linear equation 



F+[*Z± + 1)p:±zLf = (34) 

y 2/ y 2 



with the initial condition F(y) ~ y near y = 0. Clearly, the solution F(y) is a positive 
monotonically increasing function converging to a constant at infinity. The explicit solution 
is 

F(y) = ye-~M I—, — , , (35) 
where M(a, b, y) is the Kummer confluent hy per geometric function. Using the asymptotic 

■ (*)' 

small A we have 

'd+2 



expansion M(a,b,x) ~ ^x a ~ b e x for large x [15| we get F(oo) = 2T(*&)/T(&±), thus for 



B ^^-2 + Y{Af A - (36) 

For large A, repeating the argument leading to Eq. (120j) . we get for 1/A y 1 



F(y) + A ^ay ^ sin (ulny + win A + 5) . (37) 



On the other hand, for y ^> 1/A we can write 



F(y) ~ | + fT(y) , (38) 



where H(y) is a solution of the linearized equation 



H» + ( d -^ + y -) H > + d -^H = (39) 



satisfying H(oo) = B. In contrast to shrinkers, this normalization condition does not 
determine the solution uniquely since the two linearly independent solutions at infinity are 

and iJ 2 ( 2 /)~y- d exp(- 2 / 2 /4), (40) 

hence 

H(y) = BH 1 (y)+cH 2 (y), (41) 



10 



where c is an arbitrary constant. For 1/A <C y <C 1 the solutions H\ and H 2 behave as 

Hi{y)~Ciy-^wn{u\n.y + &i), % = 1,2. (42) 
Combining Eqs . (13 7 j) . f 14 1 j) . and f|42|) we get the following matching condition 

A~^a sm(u\ny + u\nA + 5) ~ BC X sin(wln?/ + Ai) + cC 2 sin (ulny + A 2 ) , (43) 
which yields 

B{A) ~CA-^sm(uj\nA + 5) } (44) 
where C and 5 are determined by a, 5, Ci, A*. 

LINEAR STABILITY OF SELF-SIMILAR SOLUTIONS 

Now, we turn our attention to the linear stability analysis of shrinkers and expanders. 
The results of this analysis are important in understanding the dynamics of blow-up. 

Shrinkers 

Substituting f(s, y) = f n {y) + w(s, y) into Eq. ffTTj) and retaining only linear terms in w, we 
get the evolution equation for linearized perturbations around the shrinker f n 

1 d 1 

w s = - (pw y ) y — cos(2/ n ) w , (45) 

which after separation of variables, w(s,y) = e~ Xs v(y), yields the eigenvalue problem 

A n v = Xv , A n = -~d y {pdy) + —^r cos(2/„) . (46) 

For each n the operator A n is self-adjoint in the Hilbert space X = L 2 ([0, 00), pdy). Both 
endpoints y = and y = 00 are of the limit-point type with admissible solutions behaving 
as v(y) ~ y for y — > and v(y) ~ y 2X for y — > 00. Note that for each n there is an eigenvalue 
A = — 1 with the associated eigenfunction v(y) = yf' n {y). The presence of this gauge mode 
is due to time translation symmetry. To see this observe that if the blow-up time is shifted 
from T to T + 2e, then 

m / ( vT j^ ) = _ " eS2//,(y) + ° {£2) - (47) 
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Since f' n {y) has (n — 1) zeroes, it follows from the Sturm oscillation theorem that for the 
n-th shrinker there are exactly (n — 1) eigenvalues below —1. We checked numerically (but 
were unable to prove analytically) that there are no eigenvalues in the interval — 1 < A < 0. 
Denoting the spectrum by {A^ n \k = 0, 1, . . . } we thus have 

< X? < ■■■ < - -K0 < < ■■■ (48) 

Concluding, the shrinker f n has exactly (n — 1) unstable modes (the gauge mode with A = — 1 
is not counted as a genuine instability). In particular, the shrinker f\ is linearly stable and 
therefore it is expected to participate in the generic dynamics of blow-up. This expectation 
will be confirmed numerically in section 6. The first few eigenvalues of the operator A n for 
several n in d = 3 are displayed in Table II. Note that the columns in this table converge to 
to limiting values, namely for each integer m we have 



hmAS m = A m . (49) 

n— >oo 



Now, we will show that X m are the eigenvalues of the point spectrum of the operator 

X d 1 

Ac = —Oy (pOy) ~ , (50) 

which is obtained from fT46l) by taking the (nonuniform) limit f n (y) — > 7r/2 as n — >■ oo. The 
potential term in ( 1501) is unbounded from below as y — > and y = is a limit-circle point, 
so for Aoo to be self-adjoint, we have to specify an additional boundary condition at y — 
(which is usually referred to as the self-adjoint extension). This is done as follows. The 
solution of the eigenvalue equation AooV = Xv which is admissible at infinity (i.e., behaving 
as v (y) ~ y 2X for y — > oo) reads 

v(y) = y 1 -^u(±-± + ^ - A, 1 + , (51) 

where U(a, b, z) is the Tricomi confluent hypergeometric function. Using the asymptotic 
expansion formula for 2; — ?- (which is valid for 1 < Re(6) < 2) jl5| 

1 [a — — 1) 1 (a) 

we get from (|5T|) 

v {y) ~ Z/ 1 ^ c °s {oolny + $(A)) as y -»• , (53) 
12 



where 

The self-adjoint extension amounts to fixing the phase $(A) modulo it. A natural choice is 
to require that the eigenvalue A = — 1 belongs to the spectrum of Aoo- This leads to the 
quantization condition 

$(A OT _i) = $(-1) + run, me Z. (55) 

As shown in Table HTl solutions of this equation, in fact, give the limit of the point spectra 
of the operators A n for n — > oo. 



n 


A n-4 


\(») 
A n-3 


X M 


\(») 
A n-1 


X (n) 
An 


\(») 

A n+1 




1 








-1 


0.51762 


1.63038 


2.69684 


2 






-53.2995 


-1 


0.48625 


1.61122 


2.68550 


3 




-6054.92 


-52.4152 


-1 


0.48271 


1.60879 


2.68380 


4 


-699295 


-5968.91 


-52.3292 


-1 


0.48237 


1.60858 


2.68363 










-1 








oo 


-688498 


-5959.55 


-52.3200 


-1 


0.48234 


1.60852 


2.68361 



TABLE II. The first few eigenvalues of the operator A n in d = 3. Numerical solutions of the 
quantization condition (155 1) are listed in the last row. 



Expanders 

The linear stability analysis of expanders proceeds along the similar lines as above. Sub- 
stituting F(a,y) = F(y) + W(a, y) into Eq.()28l) and linearizing we obtain the evolution 
equation for linearized perturbations around an expander F(y) 

W «=Ji( R W v)y ~ cos ( 2F ) W > ( 56 ) 

which after separation of variables, W(a,y) = e~ Aa V(y), leads to the eigenvalue problem 



BV = AV, B = -^-d y (Rd y ) + d —^cos(2F). (57) 

it y 
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The operator B is self-adjoint in the Hilbert space Y = L 2 ([0, oo), R dy). Both endpoints are 
of the limit-point type with admissible solutions V(y) ~ y for y — > and V(y) ~ y 2X - d e -y 2 / 4 
for y oo. The gauge mode due to time translation symmetry V(y) = yF'(y) has the 
(formal) eigenvalue A = 1, because if T — > T + 2e, then 



F(y) F ( Vl _ y 2ee _ g ) = *Xv) + eeV&l) + 0(e 2 ). 



(58) 



The gauge mode is not an eigenfunction (because it does not belong to Y), nevertheless the 
Sturm oscillation theorem still applies and implies that an expander with n zeros of F'(y) 
has exactly n eigenvalues below +1 (this was proved independently in Q|). In particular, 
monotonic expanders are linearly stable. Although this fact will be sufficient for the analysis 
of continuation beyond the generic blow-up, we wish to point out that using the "turning- 
point" method [16( one can determine sharp stability intervals for expanders. This is done 
as follows. Let F^(y) denote the expander starting with F'(0) = A. Differentiating Eq. (j2"9"|) 



with respect to A we find that 9F q^ is the zero mode of the operator B. In general, 
dF Q^ ~ B'(A) 7^ for y — > oo, however it follows from Eq. (|44|) that there is an increasing 



OA 

hence the zero mode is a genuine eigenfunction. By 16( this implies that Ak are turning 
points at which the expander Fa picks a new unstable mode. More precisely, the expander 
Fa with A G (Ak-i, Ak) has exactly (k — 1) instabilities (here Aq = by definition). 



sequence of numbers Ak (k G N) for which B'(Ak) = and then, by Eq. (1331) 

c'(A)y- d exp(-y 2 /A), (59) 



CONTINUATION BEYOND BLOW-UP 

Suppose that the solution of Eq. d^J) develops a type I singularity at time T. Then, as we 
showed above, the profile of blow-up is given by one of the shrinkers 

lim u(t, rVT - t) = fJr) , neN. (60) 

t/T 

In order to continue the solution beyond blow-up, for times t > T we glue an expander 
which matches the shrinker f„ at time T, that is we require that 

lim u{t,rVt^T) = F A {r) , F A (oo) = / n (oo) <^ B{A) = b n , (61) 

or 

\\mu(t,rVt^T) = n-F A (r), tt - F A (oo) = f n (oo) ^ B(A) = -b n , (62) 
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Note that in the case fl6T|) the solution stays continuous across blow-up (hence the degree 
does not change), while in the case (1621) the solution jumps at the origin from u(t, 0) = for 
t < T to u(t, 0) = 7r for t > T (hence the degree changes by one). In both cases we obtain 
a global weak solution which is smooth except for the time T. 

Let N(n) denote the number of roots of the equation |-B(>1)| = \b n \. It follows from the 
large A formula for expanders (jUJ) and large n formula for shrinkers (|26|) that N(n) increases 
indefinitely with n. More precisely, we find numerically (see Fig. 3) that for n > 2 



N(n) 



2n-3 for d = 3, 4, 
2n — 1 for d = 5, 6. 



(63) 



Since the shrinker f n has (n — 1) instabilities, all n > 2 blow-ups are non-generic phenomena 
of codimension (n — 1). It follows from the stability analysis of expanders that only one 
continuation is stable, namely that with A* n = min{v4 : |-B(A)| = |6 n |}. For this stable 
continuation the degree changes by one if n is odd and does not change if n is even. 
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FIG. 3. Plots of |-B(t4)| in the log- log scale. The intersections with horizontal lines \B\ = \b n \ 
determine the number of continuations beyond blowup. 
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Hereafter, we focus on the most important and interesting case n — 1 corresponding to 
the generic blow-up governed by the linearly stable shrinker f\. In this case the equation 
B(A) = bi has no roots, while the equation B(A) = —b\ has exactly one root (note that 
the existence of this root is guaranteed by the small A formula (136]) and the continuity of 
the function B(A)), hence the continuation beyond blow-up is unique, stable, and changing 
degree. As this is our main result, let us phrase it in the form of a conjecture: 

Conjecture 1. Let 3 < d < 6. Suppose that u(t,r) is a generic solution of Eq.(j9]) which 
develops a singularity at r = in a finite time T. Then, for sufficiently small r there holds 



u(t, r) ~ < 



2 +i " 



7T 



-Fa. 



for T - r 2 < t < T, 

for t = T, 
for T < t < T + r 2 . 



(64) 



where A\ is the (unique) root of equation B(A) = —b\. This is illustrated in Fig. 4. 
In the next section we present numerical evidence supporting this conjecture. 





o 1 1 1 1 1 1 1 1 1 1 1 

01 23456789 10 

y 

FIG. 4. Gluing together the stable shrinker and expander in d = 3. Here b\ = 0.573141 and 
A\ = 0.483668. 

NUMERICAL EVIDENCE 

In this section we verify the above heuristic predictions by numerical simulations. In 
order to keep track of the structure of the singularity developing on a vanishingly small 
scale, it is necessary to use an adaptive method which refines the spatio-temporal grid near 
the singularity. Our numerical method is based on the moving mesh method combined 
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with the Sundman transformation, as described in 17| . with some minor modifications and 
improvements specific to the problem at hand. This method is particularly efficient in 
computations of self-similar singularities. To implement the adaptivity in time we introduce 
a new computational time variable r defined by 



u, 



(65) 

r=0 



~r~ = g{u) , g{u) = 

ar u r t 
Under this rescaling (called the Sundman transformation) the fixed time steps in r corre- 
spond to Ati ~ (T — tj)Ar as t /* T. In this way, the time scale of the developing singularity 
is identified automatically even though the blow-up time T is unknown beforehand. To im- 
plement the adaptivity in space we introduce a new computational spatial variable £ G [0, 1] 
and define a mesh function r(£, t) which places the moving mesh points at ri(t) = r(zA£, t). 
The function r(£, t), whose role is to cluster the mesh points near the singularity, is deter- 
mined by an auxiliary moving mesh partial differential equation (MMPDE), which is solved 



simultaneously with the original PDE. We use the so called MMPDE6 



18| 



er m = -(M n )s, (66) 



with the mesh density M — \u r \ + ^/\u rr \ and the time-dependent relaxation parameter 
e(t) = 100a/ g(t) + 0.05 (this e(t), found empirically, results in a better performance than 
the customarily used constant value). 

The harmonic map heat equation ut = N(u), where N(u) is the right hand side of Eq.(j4]) 
or Eq.([9]), is now be rewritten as the system 

U = g{u) , (67) 
u T + r T u r = g(u)N(u) (68) 
er T ^ = -g{u){Mr^. (69) 

These equations are discretized using a 5-point finite difference scheme and integrated via 
the Embedded Runge-Kutta-Fehlberg (RKF45) method. 

The numerical results are presented below for d = 3 as an illustration; the behaviour of 
solutions is qualitatively the same in all dimensions 3 < d < 6. Since the dynamics of blow- 
up does not depend on the curvature of the domain, we first show simulations for Eq.([9]), and 
only at the end we show simulations of multiple blow-ups for the spherical domain Eq.(J3]). 

We begin by demonstrating the convergence to the stable shrinker. Fig. 5 depicts snap- 
shots from a typical evolution ending in a singularity. As the blow-up is approached, the 
solution is seen to converge to the profile of the stable shrinker f x . 
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According to the linearized stability analysis the deviation of the solution from the stable 
shrinker is expected to have the following form near r = for t /* T (to avoid notational 
clutter, hereafter we drop the superscript (1) on the eigenvalues and the eigenfunctions) 

oo 

u{t,r) - My) ~ ]Tc fe (T - t) A * v k {y) = Cl (T - t) Al + O ((T — t) A2 ) , (70) 
fc=i 

where y = r/y/T — t. To verify this prediction we proceed as follows. Differentiating ( ITU]) 
twice and using the normalization v[(0) = 1 we obtain 

d t [{T-t)^u r ] \ r=0 = -CL\i{T-t)*+O((T-t) > *) . (71) 

Fitting the right-hand side of this equation to the numerically computed left-hand side, 
we get the coefficient c\ and the eigenvalue Ai (see the left panel of Fig. 6). The fit gives 
Ai = 0.519, in good accord with the linearized stability analysis (see Table II). Next, in the 
right panel of Fig. 6 we show that near the blow-up time the left- and the right-hand sides 
of the expression ( 1701) (computed completely independently) do indeed agree. 

Next, we describe the continuation beyond blow-up. In order to pass through the singu- 
larity we need to modify the numerical code. First, according to (1621) we expect that at the 
blow-up time the solution is discontinuous at r = 0. This behaviour is not compatible with 
the boundary condition u(t, 0) = implemented in our code. To go around this difficulty, 
we simply rewrite Eq.Q in terms of z(t,r) = ru(t,r) and impose the boundary condition 
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FIG. 6. Left: The log-log plot of the left-hand side of expression (|7ip . The linear fit gives 
Ai = 0.519. Right: We plot the deviation of the numerical solution from the stable shrinker at 
T— t = 3.19- 1CP 6 and superimpose the first stable eigenmode c\(T — t) Xl vi(y), obtained by solving 
the eigenvalue equation (|46p . with the coefficient c\ taken from the fit in the left panel. 

z(t,0) = (which is compatible with the jump). In the case of Eq.(jH) we use a similar 
trick introducing Z(t,9) = sm(8) U(t,9) as an independent variable. Second, at some late 
stage of blow-up (say, T — t = 10 -10 ) we must switch off the Sundman transformation since 
otherwise the time step would keep decreasing down to the machine precision, effectively 
freezing the simulation and preventing it to cross the time of blow-up. To this end, we 
replace g(u) in (1651) by G(u) = g(u) + A where A 10 -10 serves as a small scale cut-off. 
When g{u) <C A, the solver loses its ability to adapt the time step appropriately and very 
quickly steps over the blow-up time. A moment afterwards, when g(u) exceeds A again, 
the Sundman transformation is turned back on and keeps tracking of the, now growing, 
time-scale of the expander. During a short time interval T — 10~ 10 < t < T + 10~ 10 when 
the time adaptation procedure is suspended, the spatio-temporal scales are unresolved and 
the numerical solution is inaccurate (the third and the fourth snapshot in Fig. 7). 

Applying this method, we continue the evolution shown in Fig. |5]past the singularity. In 
accord with Conjecture 1, almost immediately after the blow-up the numerical solution takes 
the form of the expander 7r — Fa*(u) (see Fig. 7). As written above, numerical evolution 
through a singularity necessarily involves an interval of uncontrolled behaviour due to the 
inevitable loss of resolution near the instant of blow-up. For this reason the simulation has 
limited reliability and taken alone would not provide ample evidence for the conjectured 
behaviour. It is the excellent consistency between numerics and the analytic insight, based 
on the understanding of self-similar solutions and their linear perturbations, which makes 
us feel confident that our conjecture is true. 
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FIG. 7. The same evolution as in Fig. [5] but using the modified numerical method which allows the 
solution to pass through the singularity. An additional dotted line shows the expander tt — F^* (y) . 
Notice that for \T — t\ < 10~ n the spatio-temporal resolution is lost and the numerical solution 
slightly deviates from the shrinker fi(y) (the third snapshot) and the expander n — Fa*{v) (the 
fourth snapshot). When the resolution is regained, the solution converges to the expander but 
later it moves away from it due to the interference with the far-field structure. 



Finally, let us consider the heat flow for harmonic maps between spheres U : S d — >• S d . 
As emphasized above, the curvature of the domain manifold is irrelevant in the formation 
of point singularities, hence all the above results concerning the asymptotic dynamics of 
blow-up (in particular Conjecture 1) remain valid in the case of a spherical domain. What 
makes the spherical domain interesting is a pattern of multiple blow-ups for high-degree 
initial maps. This is illustrated in Fig. 8 showing three consecutive blow-ups at the north 
pole, south pole, and again the north pole (animated simulations can be found at 3]). At 
each blow-up the degree of the map changes by one and eventually the solution comes to 
rest at the zero energy constant map. Note that, in view of the monotonicity formula (jHJ) 
and Struwe's theorem (asserting that for harmonic maps between compact manifolds the 
heat flow starting from an initial map with nonzero degree and sufficiently small energy 
must blow up in finite time), Conjecture 1 implies that the solution starting from an initial 
map of degree k must blow-up at least k times (note that the degree of the map need not 
decrease monotonically) . 20 
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FIG. 8. The solution of Eq.Q starting from the initial map Uq(6) = sin6 + 39 of degree 3. After 
three blow-ups the map becomes topologically trivial and settles down to the constant map. 

FINAL REMARKS 

As mentioned in the introduction, the global weak peaking solutions (having the form of 
a shrinker and an expander glued together at infinity) exist for many supercritical heat flow 
equations, so it is natural to ask if these equations, similarly to the harmonic map flow, enjoy 
the uniqueness of continuation beyond the generic blow-up. We are currently investigating 
this question in the following models: 

• £-equivariant harmonic map flow: The corotational ansatz (r, 0) — > (w(r),0) is 
the special [l = 1) case of a more general £-equivariant ansatz (r, (f>) — > (u(r) , xe{4>)) , 
where xe '■ S^ 1 — > S^ 1 is an eigenmap with constant energy density k = £(£+d—2)/2. 
For £-equivariant maps Eq.([9]) changes to 

^ = ^T(^V) r -^sin(2.). (72) 

All the qualitative results concerning existence of shrinkers and expanders and their 
linear stability obtained above for £ = 1 trivially carry over to £ > 1 provided that 
3 < d < 2£ + 2y/I + 2, however the quantitative characteristics of self-similar solutions 
(in particular those which imply the uniqueness of gluing an expander to the stable 
shrinker) remain to be checked. 
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Yang-Mills heat flow: It is well-known that there are close parallels between the 



20| | . For the spherically symmetric 



harmonic map and the Yang-Mills heat flows 

magnetic Yang-Mills potential h(t,r) in d > 3 dimensions the analogue of Eq.fJH]) 
reads 

h t = -jL (r^ 3 /.,,), - d —^ L h{h- l)(h - 2) . (73) 



Using a similar shooting technique as in {3] one can easily show that for 5 < d < 9 
there are infinitely many shrinkers h(t,r) = <f) n (y). One novel feature, in comparison 
with the harmonic map flow, is that the first (stable) shrinker is known explicitly 21] : 



l(v) 



h^d - 12 - (d + 2)V2d - 4) 



x/^2 



6 + cm/ 2 ' 2 V " ' ■ _V2 

This may be helpful in proving the Yang-Mills analogue of Conjecture 1. 

Semilinear heat equation: The equation 



(74) 



u t = Am + \u\ p u 



(75) 



for d > 3 and supercritical powers 



d+2 
d-2 



< p < p :- 



1 + 



d-10 



for 3 < d < 10, 
for d > 11, 



(76) 



has self-similar solutions (shrinkers and expanders) 22|-|24j| which give rise to global 



peaking solutions similar to the ones described in section 5, however all these solutions 
are unstable (cf. also 25, [26]). It seems interesting to see if a kind of analogue of 
Conjecture 1 holds for (codimension-one) threshold solutions. 

In this paper we restricted our analysis to dimensions 3 < d < 6. We wish to emphasize 
that this is not a technical restriction. For d > 7 the shrinkers disappear and consequently 
the blow-up changes character from type I to type II j^ . 
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